Well-established nucleon resonances revisited by double-polarization measurements 
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The first measurement is reported of the double-polarization observable G in photoproduction of 
neutral pions off protons, covering the photon energy range from 620 to 1120 MeV and the full solid 
angle. G describes the correlation between the photon polarization plane and the scattering plane for 
protons polarized along the direction of the incoming photon. The observable is highly sensitive to 
contributions from baryon resonances. The new results are compared to the predictions from SAID, 
MAID, and BnGa partial wave analyses. In spite of the long-lasting efforts to understand •yp — > pir° 
as the simplest photoproduction reaction, surprisingly large differences between the new data and 
the latest predictions are observed which are traced to different contributions of the A(1535) with 
spin-parity J p — 1/2 - and A(1520) with J p — 3/2~. In the third resonance region, where JV(1680) 
with J p — 5/2 + production dominates, the new data are reasonably close to the predictions. 

PACS numbers: 14.20 



Symmetry arguments have led Cell-Mann and Zwcig to 
introduce the concept of quarks as constituents of mesons 
and baryons 0, Q. Later it was realized that quarks 
carry a new kind of charge, color, and are the sources 
of gluons transmitting the binding forces. A new theory 
evolved, quantum chromodynamics [H, which proved to 
be very successful for large momenta where the interac- 
tion becomes weaker. In the region of meson and baryon 
resonances, QCD resists any perturbative approach, and 
numerical calculations on a space-time lattice are nec- 
essary. It has been only recently that the full baryon 
spectrum including physical quantum numbers was pre- 
sented Q , even though unrealistically large quark masses 
had to be used. Qualitatively, these lattice results con- 
firm the classic quark model calculations Q , whereas the 
agreement with a fully relativistic model is less convinc- 
ing Yet the experimental mass spectrum exhibits 
some remarkable differences to these calculations. The 
AT (1440) resonance with spin-parity J p = l/2 + has a 
mass of 1440 MeV and falls below its spin-parity partner 
A(1535) with J p = 1/2" while QCD on the lattice and 
quark model calculations predict the reverse. Also, above 
a mass of 1800 MeV, the number of observed resonances 



falls short of the number of predicted resonances. 

Baryon resonances have very short life times, hence 
many resonances with different spin-parities overlap; 
their properties have to be unfolded from data in partial 
wave analyses. Most partial wave analyses are done us- 
ing ttN elastic and charge exchange scattering, but pho- 
toproduction experiments offer the opportunity of ana- 
lyzing resonances which were not accessible up to now. 
However, the number of independent observables neces- 
sary for a unique partial wave solution increases. Below 
the two-pion production threshold at E 1 = 310 MeV, the 
Watson theorem [?| relates the irN phases to those in 
photoproduction of single pions. In this case the mea- 
surement of the differential cross section da and of the 
photon beam asymmetry S is sufficient to determine pre- 
cisely the s- and p-wavc multipolcs (Z ff = 0, 1), as well as 
the small ratio of electric quadrupole (E2) and magnetic 
dipole (Ml) amplitudes in the N A(1232) transition 
d, [§]. Above the two-pion threshold, but still in the 
mass range where nucleon resonances carry only one unit 
of intrinsic orbital angular momentum, five experiments 
are sufficient to construct the spin-helicity structure of 
the photoproduction amplitude 13, 11 1. In general, the 
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Figure 1. (Color online) a) The central part of the detector system. The CsI(Tl) crystals (blue and green) are read out via 
wavelength shifters and photodiodes or photomultipliers, BaF2 crystals (yellow) in forward direction with photomultipliers. 
b) The missing mass distribution, and (c) the 7r° angular distribution for reaction (1) for an incident photon energy of _E 7 = 
1000±25 MeV; butanol (□), hydrogen (A), carbon (o), and the sum of hydrogen and carbon data (*). From these distributions, 
the dilution factor (d) is determined. 



observation of eight polarization observablcs with high 
precision is required to construct the electric (E) and 
magnetic (M) multipoles governing the process fl2i. Il3|. 

In this letter we report on a first measurement of the 
double-polarization observable G in photoproduction of 
neutral pions off protons covering the second and third 
resonance region. Linearly polarized photons we re p ro- 
duced by scattering of a 3.2 GeV electron beam [14| off 
a diamond crystal [Tij . whereby maximal polarization 
of 65% at 950 MeV and 59% at 1150 MeV for a sec- 
ond data set were reached. The photons then hit a bu- 
tanol (C4H40O) target with longitudinally polarized pro- 
tons [l6j, with a mean proton polarization of about 75%. 
The butanol target was replaced by a hydrogen or car- 
bon target for background studies and for normalization. 
The incoming photons may produce a 7r° in the reaction 
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(1) 



For linearly polarized photons (with polarization p 7 ) 
and protons (with polarization pt), the number of events 
N at the polar angle 9^ due to reaction ([1) as a function 
of the azimuthal angle 0^ can be written in the form 



= 1 - p 7 £ B cos(2^) + p 7 p T G B sin(2^) (2) 



where N is given by averaging N(^ n ,0 n ) over <fr„. The 
beam asymmetry Eb and the observable Gb for the bu- 
tanol target are related to the corresponding quantities 
for scattering off free (/) protons S, G and bound (b) 
nucleons by 



E R = 



N[ S + V b E b 



Gb 



Nn + N't 



G. (3) 



The number of events in the denominator is the sum 
of the number of events where scattering took place off 
protons or a bound nucleon, N n = Nq + Nq. There is 
no contribution of bound nucleons to Gb since carbon 
and oxygen nuclei carry no polarization. The quantity 
D = iVg /(Nq + JVq ) is called dilution factor and depends 
on the reaction and on the kincmatical cuts performed. 



It is also possible to determine G without using the di- 
lution factor by reversing the target polarization and re- 
garding the difference N + ((f>„, 9^) — N~(4' 7T , 6 W ), given by 
2Nq p 7 pxGsin(2</) T ). A similar equation can be written 
for a change of the photon polarization plane by ±7r/4 
171. 



Fig. Q]a, shows the calorimeter setup. Neutral pions 
are reconstructed from their decay into two photons and 
a measurement of their energy and direction in CsI(Tl) 
and BaF2 crystals enclosing the target hermetically. The 
proton direction is determined, assuming that it origi- 
nated in the target center, from its hit in the surrounding 
three-layer scintillation fiber detector [l8[ or scintillation 
detectors in front of the crystals in forward direction, 
and its hit in the Csl or BaF2 crystals. In comparison 
to (l9| . two major instrumental changes have been in- 
troduced. The forward opening of the calorimeter is now 
covered by 90 additional CsI(Tl) crystals with photomul- 
tiplier readout in the main calorimeter and 216 BaF2 
crystals, 2.1m downstream of the target. The tagger is 
now mounted horizontally. Beam electrons not produc- 
ing brcmsstrahlung arc now deflected to the left and are 
stopped in a beam dump well behind the detector system. 

In the first step of the reconstruction, events due to re- 
action (1) are selected using a series of kinematical cuts. 
Two classes of events are retained, events with three hits 
in the calorimeters and events with two hits. Both hits 
of a two-hit event are assumed to be photons; their in- 
variant mass is calculated and required to fall into a 3<r 
(±25 MeV) window centered at the 7r° mass. The proton 
is assumed to be stopped in some material, but due to 
energy and momentum conservation, the proton momen- 
tum can be reconstructed. If there is a third hit, one pair 
of hits is assumed to be photons, and it is checked that 
the pair forms a 7r°, using again a ±3<r window. The 
third hit has to match the direction of the proton mo- 
mentum reconstructed as missing momentum. Matching 
is defined by a cone, ±7° in the azimuthal and polar an- 
gles. Finally, a time coincidence is required between the 
tagger hit and hits in one of the scintillation detectors or 
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The beam asymmetry E 
23 



Figure 2. (Color online) a) A typical (j)^ -distribution with a fit using eq. (2). b,c) 

of cos#,r for E-y — 800 MeV and for E-, — 1100 MeV. Black dots show our data, red squares GRAAL data 
represent predictions from different partial wave analyses. Solid (black) curve: BnGa [24||; dashed (red): SAID 
(black): BnGa with E + and E 2 - amplitudes from SAID; dotted (blue): MAID [2j|; dashed-dotted (black): BnGa with E + 
and E 2 - amplitudes from MAID. Gray area shows the systematic error due to interactions on nuclei and uncertainty in the 
photon polarization. 



as a function 
. The curves 
; long-dashed 



forward calorimeter crystals. Additionally all events are 
removed, if the calculated beam photon energy (under 
the assumption that the reaction occurred on a proton) 
falls below the tagging threshold. With these cuts the 
reconstruction efficiency for reaction ([1} on free protons 
in the butanol target exceeds 60%. 

Fig. [T}d shows the resulting distribution of missing 
masses for events passing the selection criteria, with bu- 
tanol, carbon, or hydrogen as target material. By proper 
scaling, the sum of distributions from hydrogen and car- 
bon reproduces very well the distribution obtained with 
the butanol target. The missing mass distributions show 
a peak at the proton mass and a tail due to events where 
an additional meson was produced but escaped detec- 
tion. A final cut is now applied on the missing mass at 
nip ± 50 MeV. The resulting cos 8„ distribution of these 
events is exhibited in fig. [TJ;. Again, the sum of distribu- 
tions from hydrogen and carbon match the distribution 
obtained with the butanol target; there is no new scaling 
factor involved. 

For the double-polarization observable G the correc- 
tion factor for reactions on nuclei is the dilution factor, 
displayed in fig. QJL Due to the kinematic cuts, the di- 
lution factor depends on both the photon energy E 7 and 
the 7T° scattering angle 9 n . For a backward 7r°, the proton 
momentum is high and therefore detected in the crystals 
with a high efficiency. 

A typical distribution of the measured (/^-dependent 
data yields with a fit using eq. (2) is shown in fig. [2ji. 
Obviously, there is a large cos 20^ contribution but the 
data are not symmetric with respect to <p v = 0: there 
is a significant sin 20^ contribution from which the de- 
sired observable G is deduced. Fig. [5b, c show the beam 
asymmetry E extracted from this and other similar plots 
for two energy bins. The systematic error due to reac- 
tions on nuclei are shown as gray bars. The new E data 
are in excellent agreement with earlier measurements of 
GRAAL 0. 

The observable G is determined from eq. ©, as well 




Figure 3. (Color online) The polarization observable G as a 
function of cos 6„ from _E 7 = 800 MeV up to E 1 = 1100 MeV. 
Systematic errors are shown in gray bars. Curves: see fig. [2] 



as by reversing the target polarization, or changing the 
photon polarization plane by ±7r/4. The results are con- 
sistent, their spread is used to determine the systematic 
errors. These include the systematic uncertainty in the 
dilution factor, errors due to background in cos^ > 0.9, 
the relative uncertainty in the target polarization of ±2%, 
and the relative uncertainty in the photon polarization 
of ±5% [2l|. Fig. [3] shows the angular distribution of 



the polarization observable G for selected photon energy 
ranges. The error bars give the statistical errors, the 
systematic errors are shown as gray bars. 

The new data on E and G are compared to the predic- 
tions from the SAID (SN11) 0, MAID [H, and BnGa 
[24] ] partial wave analyses. The results are surprising. 

Already at rather low energies, in the region of the 
four-star resonances V(1440), V(1535) and V(1520), 
very significant differences in the predictions can be ob- 
served (see figs. [2] and [3]). At 800 and 900 MeV, the 
SAID prediction is incompatible with the data on E and 
G, MAID disagrees with the data on E and on G at 
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Figure 4. The double-polarization observable G as a function 
of energy for two selected bins in 6^. Curves: see fig. [2] For 
comparison the positions of several resonances are marked. 



900 McV. At 1000 and 1100 McV, BnGa, MAID and 
SAID give a reasonable description, even though MAID 
is higher at 1000 McV. These conclusions are supported 
by fig. 0] where G is shown as a function of E 1 for two 
scattering angles and compared to PWA predictions. At 
low energies, SAID deviates strongly from the data; in 
the 900 — 1000 McV range, the consistency between data 
and MAID or SAID predictions is poor. Above, none of 
the partial wave analyses reproduces the data. 

We have studied possible reasons for the discrepant re- 
sults. In this discussion, the two electric multipoles E + 
and E 2 - play a leading role. If we replace the BnGa E + 
and E 2 - multipoles by the corresponding SAID multi- 
poles, we reproduce the SAID results for £ and G over a 
wide energy range, and the same observation holds true 
for MAID. The reason for the deviation between data and 
MAID and SAID is obviously traced to these two multi- 
poles. All other multipoles can be exchanged without a 
large impact. 

Resonances with no angular momentum between pro- 
ton and 7T° (l n = 0) contribute to the E + multipolc 
which is sensitive to the iV(1535), A(1620) and V(1650) 
resonances. The resonances A r (1520) and A(1700) with 
1-^ = 2 and total spin J = 3/2 contribute to the E 2 - 
multipolc. Indeed, the three partial wave analyses BnGa, 
MAID, and SAID give significantly different hclicity am- 
plitudes for these resonances, in particular for the Ai/ 2 
amplitude of A(1535) and for A 3/2 of A(1520) (see Ta- 
ble©. 

In the medium energy range, in the third resonance 
region, the jp — > pir° cross section is dominated by the 
J p = 5/2 + V(1680) resonance. In this regime, there is 





Af(1535)l/2~ 


iV(1650)l/2~ 


A(1620)l/2~ 




Solution 


A 1 ' 2 


A 1 ' 2 


A 1 ' 2 




BG2011-02 


105 ± 10 


33 ± 7 


52 ± 5 




MAID-2007 


66 


33 


66 




SAID-2011 


99 ± 2 


65 ± 25 


64 ± 2 






JV(1520)3/2~ 


A(1700)3/2~ 


Solution 


A 1 ' 2 


A 3 ' 2 


A 1 ' 2 


A 3 ' 2 


BG2011-02 


-22 ± 4 


131 ± 10 


160 ± 20 


165 ± 25 


MAID-2007 


-27 


161 


226 


210 


SAID-2011 


-16 ± 2 


156 ± 2 


109 ± 4 


84 ±2 



Table I. Helicity amplitudes of low-lying negative-parity N 
and A resonances contributing to the E + and E 2 - multipoles 
(in GeV-i 1(T 3 ). 



good agreement between the three predictions. The prop- 
erties of this resonance are obviously well defined: For 
the leading A 3 / 2 helicity amplitude, BnGa finds 135 ±6, 
MAID 134, and SAID 141 ± 3 (in GcV^ 10~ 3 ). 

In the high energy region, increasing differences be- 
tween the predictions can be observed for some angles. 
This is the so-called fourth resonance region where many 
resonances exist with questionable evidence for their exis- 
tence and, if they exist at all, with at least poorly defined 
properties. New data on G in this energy region will help 
to disentangle the spectrum of nucleon resonances [2ll ] . 

Summarizing, we have reported the first measurement 
of the double-polarization observable G in the reaction 
jp — > pir° in a wide range of energies and covering the full 
solid angle. The new data on G resolve the discrepant re- 
sults on the helicity amplitudes of low-lying four-star nu- 
cleon resonances obtained from BnGa, MAID, and SAID 
and are an important step towards a complete data base 
which will define unambiguously the nucleon excitation 
spectrum. 

We thank the technical staff of ELSA and the parti- 
cipating institutions for their invaluable contributions to 
the success of the experiment. We acknowledge support 
from the Deutsche Forschungsgemeinschaft (SFB/TR16) 
and Schweizerischer National] 'onds. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 
[10] 

[11] 
[12] 

[13] 



Gell-Mann, Phys. Rev. 125, 1067 (1962). 
Zweig, CERN-TH-401 (1964). 

Fritzsch, P. Minkowski, Nuovo Cim. A30, 393 (1975). 
G. Edwards et al, Phys. Rev. D84, 074508 (2011). 
Isgur, G. Karl, Phys. Rev. D18, 4187 (1978). 
Loring et al, Eur. Phys. J. A10, 395 (2001). 
M. Watson, Phys. Rev. 182, 228 (1954). 
Beck et al, Phys. Rev. Lett. 78, 606 (1997). 
Blanpied et al., Phys. Rev. Lett. 79, 4337 (1997). 
S. Omelaenko, Yad. Fiz. 34, 730 (1981). 
Y. Wunderlich, Diploma Thesis, Univ. Bonn (2012). 
W. -T. Chiang, F. Tabakin, Phys. Rev. C55, 2054 (1997). 
A. M. Sandorfi et al., J. Phys. G 38, 053001 (2011). 



5 



[14] W. Hillert, Eur. Phys. J. A28S1, 139 (2006). 
[15] D. Eisner et ai, Eur. Phys. J. A39, 373 (2009). 
[16] H. Dutz et al, Phys. Rev. Lett. 93, 032003 (2004). 
[17] A. Thiel, PhD Thesis, Universitat Bonn (2012). 
[18] G. Suft et al, NIM A538, 416 (2005). 
[19] V. Crede et ai, Phys. Rev. C84, 055203 (2011). 



[20] O. Bartalini et ai, Eur. Phys. J. A26 399 (2005). 
[21] A. Thiel, H. Eberhardt et al., publication in preparation. 
[22] R. Workman et al. Phys. Rev. C85, 025201 (2012). 
[23] D. Drechsel et al, Nucl. Phys. A645, 145 (1999). 
[24] A. V. Anisovich et ai, Eur. Phys. J. A48, 15 (2012). 



